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The degree of ionisation of a two-dimensional electron-hole plasma is calculated in the low-density 
(Boltzmann) limit. The electron-hole interaction is considered for all states: optically active and 
^ , inactive, bound and unbound. The theory is applied to exciton/free-carrier plasma in GaN-based 

^ ■ quantum wells at room temperature. 

' Room-temperature operation of blue quantum- well lasers, based on GaN Q or ZnSe [^, is no longer a novelty. 
However, the nature of lasing in wide-gap semiconductors and quantum wells has not yet been completely understood, 
^sj The large exciton binding energy in wide-gap semiconductors and quantum wells favours excitonic gain processes for 
which no satisfactory theoretical treatment exists . A knowledge of the balance between excitons and free carriers 
(-H is crucial in determining the dominant gain process. It is known that a naive application of the Mott criterion for the 
^ metal-insulator transition as well as the use of a single-bound-state law of mass action are insufficient, as screening of 
^ \ excitons by the electron- hole plasma and strong scattering of particles within the plasma both play a crucial role . 

For a consistent description of the electron-hole plasma at a temperature which is higher then the exciton binding 
energy, bound states and unbound scattering states should be treated on the same footing. In what follows we present 
i such a consistent treatment for the purely two-dimensional (2D) case in the low-density (Boltzmann) limit. 
c/2 \ Following an approach applied in 3D to nuclear matter [gj, an ionic plasma Q, and the electron- hole system in 
excited semiconductors [Q, we divide the total electron (hole) density between two terms: 
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The first term nJJ is the density of uncorrelated quasiparticles with renormalized energies. Only this term should 
be taken for the screening radius calculation Q. All correlation effects both in the bound and continuum states are 
^ ' incorporated into the second term n™''^ which is called the correlated density. The lower index in Eq. (Q) is a species 
index, a = e for electrons and a = h for holes. 
I ' In the low-density limit there is no need to go beyond two-particle correlations. This allows us to separate clearly 
J>. , the role of the inter-particle Coulomb interaction from the phase-space filling effects. It is tempting to relate n^°^^ 
and by a simple law of mass action with the single exciton bound state energy reduced by screening [^|j7|. The 
main shortcoming of this approach is a disregard of the strong scattering of unbound carriers. A complete account 
' of scattering states as well as all (optically active and inactive) bound states requires the calculation of a two-body 
partition function which involves summation over all two-particle states. In the low-density (non-degenerate) limit, 
for which there is no Pauli blocking, a 2D analogue of the modified mass action law reads 
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where p^ah — "niami,/ {nia + rrib) is the reduced effective mass, and Zab is the two-body interaction part of the partition 
function. Note that due to charge-neutrality the total electron-hole density rte = rt/i = n is independent of species, 
Q ' whereas ^ n° and n^"'"'' ^ n^"'"'' if the electron and hole have different masses. 
O ' The electron-hole part of the partition function which exhibits bound states (excitons) is given by 
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where f3 — 1/ {ksT), mfi is the projection of the angular momentum onto the axis normal to the plane of 2D motion 
(m = 0,±1,±2, . . .), h'^k^/2fj,eh is the energy of the relative motion of the unbound (scattered) electron and hole, 
k is the absolute value of the relative motion momentum, d,n{k) are the 2D scattering phase shifts introduced in 
the standard way £^,1/ are the bound-state energies (index i> enumerates bound states with given m), and the 
double sum in the first term ranges only over bound states. Equation is the 2D analogue of the Beth-Uhlenbeck 
formula |^ , and it is derived in the same fashion as in the 3D case . The scattering term in the right-hand side 
of Eq. (^ gives the contribution to Z^h of the electron-hole attraction in the continuum part of the energy spectrum. 
The electron-electron and hole-hole parts of the partition function Zee and Zhh contain the scattering term only. 
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Equations provide a consistent description of the ionisation degree, defined as a = n^/ {'nP^+n'^'^^) . Technically 
the most difficult problem is to calculate the binding energies and scattering phase shifts in a screened Coulomb 
potential. We use for this purpose the variable-phase method [|l2| known from scattering theory. In this method 
the scattering phase shift and the function defining bound-state energies can be obtained as a large distance limit of 
the phase function, which satisfies the first-order, nonlinear Riccati equation originating from the radial Schrodinger 
equation. 

In this paper we model the screened Coulomb interaction in a 2D plasma by the well-known Thomas-Fermi expres- 
sion for a statically screened Coulomb potential 

Vs{p) = T- r ^-^dq = T- (- - Jg.[Hofep) - Yoiqsp)]] , (4) 

where Qs is the 2D screening wavenumber (which depends on temperature and carrier density), e is the static 
dielectric constant of the semiconductor, Jo{x), Yo(x), and Ho(x) are the Bessel functions of the first and of the 
second kind and the Struve function. The upper sign in Eq. is for electron-hole attraction, the lower sign is for 
electron-electron or hole-hole repulsion. Being the long-wavelength static limit of the random phase approximation 
for a purely 2D case, Eq. is the simplest model for the screened Coulomb potential in 2D. Nevertheless, this 
expression reflects the fact that the statically screened potential in 2D decreases at large distances slower than in 
the 3D case. Despite numerous realistic corrections, Eq. remains the most widely used approximation for the 
2D screening, especially for the screened exciton problem. The variable-phase method application to scattering and 
bound states in the screened Coulomb potential (Q) is described in detail in our recent paper [p^ . The method is 
especially effective for calculation of shallow-state binding energies and low-energy scattering phase shifts. Applying 
the variable phase method together with a 2D analogue of the Levinson's theorem we have found that with 
decreasing Qs several bound states with different angular momenta appear simultaneously at certain integer values 
of l/(gsa*), where a* is the effective (excitonic) Bohr radius. This degeneracy is different from the well-known 
degeneracy of the unscreened 2D exciton states. In the low-density limit {qsO.* — > 0) the number of bound states 
oscillates around l/{qsa*) with the period and amplitude of oscillations proportional to 1/^/q^a*. Then, using 
the expression for the Thomas-Fermi 2D screening wavenumber qs for a two-component non-degenerate electron-hole 
plasma g^a* — An (Ry* /kBT) (n^a*^ -|- n°a*^), Ry* being the excitonic Rydberg, one can find the ratio of the 
free-carrier density to the total density. For a model semiconductor with me — mh, this ratio is equal to 2/3 in the 
low-density limit. This result is different from the single-bound-state law of the mass action, which gives a complete 
ionisation of bound states (a — > 1) in the same limit. 

Figure 1 shows the results from the calculation of the electron- hole part of the partition function, Z^h, which 
contains both the bound state sum and the scattering phase shift integral. In this figure Z^h is plotted as a function 
of the inverse screening wave number l/qs measured in units of the effective Bohr radius a*. The temperature is 
given in units of the bulk excitonic Rydberg Ry* . 
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FIG. 1. The electron-hole part of the partition function, Zeh versus the screening length l/qs for two values of fcsT/Ry*. 
Solid lines show the bound state contributions Zbound only. Dashed lines correspond to the total partition function with 
scattering states included. 



To emphasize the role of scattering we show on the same plot the bound-state sum, Zbound — J2m v 6xp(— i,). 
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which exhibits jumps whenever new bound states appear. These jumps become higher with increasing screening 
length 1 1 qs since several bound states appear simultaneously [T3| . A proper account of scattering eliminates these 
unphysical jumps. 

We calculated the density dependence of the degree of ionisation a in the 2D low-density (Boltzmann) limit 
for fcsT = 1 Ry*, which roughly corresponds to GaN at room temperature. Over a wide range of pair densities, 
0.01 < na*^ < 0.1, a is almost independent of density, and it increases outside this range. The minimal value of the 
degree of ionisation, amin ~ 0.34, corresponds to na*^ « 0.04 (n « 5 x lO-'^^cm"^ for GaN). Even at the relatively 
high densities of 10^^ cm~^ we find that there is a single bound state having a binding energy of the order of /cbT 
which is available to participate in excitonic lasing. 

We also calculated the second virial coefficient of the dilute 2D electron-hole plasma in GaAs and GaN-based 
quantum wells at room temperature. These calculations show a striking difference between a strongly correlated 
exciton/free-carrier plasma in GaN and a nearly ideal electron- hole gas in GaAs. 
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